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\ We describe black holes in d + 3 dimensions, whose thermodynamic properties correspond to 

those of a scale invariant non-relativistic d + 1 dimensional quantum system with dynamical 
exponent z = 2. The gravitational model involves a massive abelian vector field and a scalar 
field, in addition to the metric. The energy per particle in the dual theory is \fi\ d/(d+2) at 



o 



1} , any temperature, while the ratio of shear viscosity to entropy density is H/Att. 

(N 

(N 1 I. INTRODUCTION 

Scale invariance is an important concept in physics, with fundamental applications in diverse areas 
including critical phenomena and high energy physics. One can think about scale invariance as a 
symmetry under simultaneous rescaling of both space and time x l — > \x\ t — ► X z t, where A > 
is the scaling parameter, and z is the dynamical exponent. Relativistic scale invariant systems 
such as blackbody radiation have z = 1, while many non-relativistic scale-invariant systems have 
' z = 2. An interacting non-relativistic quantum system which has recently attracted much interest 

is given by a collection of cold fermionic atoms in the unitary limit In this limit the scattering 
length of atoms formally diverges whereas the mean particle distance is considered (arbitrary) large 
compared to the range of the interaction. The only relevant length scales are then given by the 
00 \ mean particle distance and the thermal wavelength, which are determined through the fermion 

density or chemical potential and the temperature. All quantities can therefore be expressed by 
their analogues in the non-interacting theory times a universal function only depending on the 
dimensionless ratio of chemical potential and temperature [d]. Consequently, the system at zero 
temperature and density should be described in terms of a non-relativistic scale invariant theory 



with z = 2 

The methods of the gauge-gravity duality [J] have been extremely helpful for understanding 
strongly interacting scale invariant relativistic models, and it is reasonable to ask whether the same 
tools can be useful in understanding strongly interacting scale invariant non-relativistic models. 
The key idea lies in realizing the relevant symmetry group geometrically, rather than directly on 
the Hilbert space of the quantum theory. The algebra of the non-relativistic conformal group 
(the Schrodinger group) can be obtained from the relativistic group by reducing along the light- 
cone and identifying the corresponding momentum with the particle number operator 0]. A 
spacetime realizing this symmetry as its isometry group was recently identified in 0] for z = 2, 
and independently in [6] for arbitrary z: 



ds 2 = -r 2z du 2 - 2r 2 dudv + r 2 dx i dx i H T . (1) 
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Here r is the radial coordinate in the standard gauge-gravity duality description with the boundary 
at r^oo, u is the boundary time coordinate, x % are d-dimensional spatial coordinates, and v is an 
additional coordinate whose conserved momentum corresponds to the particle number operator. 
For d = 2 and z = 2, this spacetime and its non-extremal generalizations were embedded in type 
IIB string theory 0, @, 0| and the equation of state in the dual field theory together with its shear 
viscosity were determined. In addition, the authors of 0, E| found a five-dimensional gravitational 
model which admits a non-extremal version of ([I]) as a solution, corresponding to a 2+1 dimensional 
scale-invariant field theory at non-zero temperature and density. 1 

Given that most experiments on cold atomic gases are performed in d = 3 spatial dimensions, 
it would be interesting to find a simple gravitational model, corresponding to a 3+1 dimensional 
scale invariant quantum system at non-zero temperature and density. Such six-dimensional grav- 
itational action would be a natural starting point for phenomenological holographic models of 
non-relativistic systems, similar to the AdS/QCD program [l3 |. It is the purpose of this note 
to present black hole spacetimes whose thermodynamic properties correspond to those of a scale- 
invariant non-relativistic quantum system in arbitrary spatial dimension d ^ 2, including the 
physically interesting case d = 3. Using the standard dictionary of the gauge-gravity duality, 
we calculate the thermodynamic potentials and the shear viscosity of the dual d+1 dimensional 
non-relativistic field theory. Our notations closely follow the presentation of Herzog et al 0]. 



II. THERMODYNAMICS OF A SCALE-INVARIANT SYSTEM 



Before presenting black hole spacetimes whose thermodynamics corresponds to that of a scale 
invariant system, let us briefly discuss what we expect thermodynamic potentials to look like. 
We are interested in a translation invariant quantum field theory in d+1 spacetime dimensions 
in the thermodynamic limit in the grand canonical ensemble given by the density operator p = 
Z^ 1 exp[— (H— pN)/T] where Z = tr exp[— (H— [iN)/T] is the partition function, T is temperature, 
H is the Hamiltonian, fi is the chemical potential, and N is the conserved particle number operator 
which commutes with H. Scale invariance means that there is a time-independent operator D whose 
commutation relations within the Schrodinger algebra include 2 

[D,H]=izH, [D,N]=i(2-z)N. (2) 

The Hamiltonian density Ti. has scaling dimension d+z, i.e. [D,TC(0)] = i(d+z)H.(0). Taking 
the expectation value of this equation in the grand canonical ensemble, using the commutation 
relations ([2]) and cyclicity of the trace gives 



1 Other recent work on holographic realization of non-relativistic theories includes [Tol . [ill . . 
The second relation in (J2j) holds in systems where particle number current coincides with momentum density. From 
a geometric perspective, this commutation relation emerges from the isometries of the spacetime |T} 6] . We thank 
J. McGreevy for a discussion on these points. 



3 



where e = tr(pH) is the equilibrium energy density. The solution is e = T~ g(T~^~ fi) where g 
is an arbitrary function. Once the energy density is known, pressure p can be determined from 
the relation e + p = Ts + un, where s = dp/dT is entropy density, and n = dp/dfj, is the number 
density. Demanding that p— >0 as e— >0, one finds for z = 2 

e(T,ti)=T*Pg(n/T), 2 e(T, a) = dp(T, a) . (3) 

We will now describe black holes in d+3 dimensions whose thermodynamics has exactly the scaling 
form of Eq. ([3]), with a particular function g(fi/T). 



III. GRAVITATIONAL MODEL 

We propose the following action of gravity coupled to a massive abelian vector field plus a scalar: 

5 = i6^^ J dd+3x V=g(R- " \^f^ - ^-A^ - v(<j>)) , (4) 

where Gd+3 is the (d+3)-dimensional Newton's constant, the scalar potential is given by V(cp) = 
(A+A')e a< ^ + (A-A')e^, and the coefficients are 

A = -l(d+l)(d+2) , A' = ^(d+2)(d+3), m 2 = 2{d+2) , a = (d+2)b = 2^- . 

For = 0, one recovers the model of refs. 0, 0], where the spacetime (pQ) emerges as a solution 
to Einstein equations with a non-trivial profile of A^. For 0^0, our model admits the following 
planar black hole solutions which are asymptotic to ([I]) 



ds = r 2 h d +! 
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(/-I) 

4(fc-l) 
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/ r V - (l+f)dudv + -^-dw 2 ) + h—i lr 2 dx l dx l + I ,(5a) 

A = ^-r 2 du + l^dv , (5b) 
2h h 

4>=--lnh, (5c) 

where h(r) = 1 + /3 2 Fq +2 /r d and /(r) = 1 — rQ +2 /r d+2 . Here (3 is an arbitrary parameter, and the 
horizon is at r = tq. It is helpful to introduce coordinates t = (3v + uj (2(3) and y = —j3v + uf (2/3), 
so that the solution © becomes 

d i dv 2 

ds 2 = r 2 h~^m (-fdt 2 + dy 2 - r 2 f(3 2 {dt+dy) 2 + hdx 1 dx l ) + -jj , (6a) 

A{r) = r ^-{fdt + dy). (6b) 

For d=2 these solutions were discussed in ref. %. Note that Eq. © with /(r) = 1, /i(r) = l + Q. d /r d 
is also a solution, with f2 a free parameter. We will now show that thermodynamics of black holes 
([5]), (|6|) matches thermodynamics of d+1 dimensional scale- invariant field theories with 2 = 2. 
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IV. BLACK HOLE THERMODYNAMICS 

The Bekenstein-Hawking entropy of the black hole is proportional to the area of the horizon, 

S=j^—r d +1 pAvV, (7) 

where we assume that v has a finite size Av, and V = Ax 1 . . . Ax d is the (infinite) volume along the 
x % directions which we identify with the spatial volume of the dual field theory. The temperature 
of the black hole can be found from the surface gravity k, 

« 2 = ~(v'*n(v^) ) (8) 



where £ is a Killing vector field which is null at the horizon. Following [8|, we choose £* = 1//3, or 

£ = - + —- (9) 

? 8u 2/5 2 8v ' W 

This gives the black hole temperature T = k/2tt equal to 

which we identify as the temperature of the dual field theory. From @ we can also infer the value 
of the chemical potential. With the generator of translations in the u-direction corresponding 
to the Hamiltonian and the generator of translations in the f-direction corresponding to particle 
number, we expect the chemical potential to be 0] 

as the conjugate Lagrange multiplier to the particle number (this is analogous to what happens in 
thermodynamics of Kerr black holes). With this identification, the chemical potential is negative. 

We identify the partition function as Z = e~ SE where Se is the the on-shell Euclidean action, 
obtained after taking t — > it. The on-shell action is divergent at large r, and requires regularization. 
To make the action finite, we add counterterms at the boundary r— >oo, 

S^S + —^ [d d+2 x^h(K -CQ + ...) , (12) 



where K is the standard Gibbons-Hawking boundary term, cq is the boundary cosmological con- 
stant, and dots denote "matter" counterterms which are polynomial in the fields (ft, and their 
derivatives. The on-shell action evaluated on the solution ((5J can be made finite by choosing 
cq = d+1, together with the appropriate coefficients of the matter counterterms, giving 3 



1 

l6irG d +3 



S E = —^ZF< l dd+2x r +2 • ( 13 ) 



3 Alternatively, the Euclidean on-shell action can be regularized by subtracting the action of the zero-temperature 
solution ([5]) with f(r) — 1 and h(r) = 1 + /3 2 rg +2 /r d , and demanding that the volumes at large fixed r coincide in 
the two geometries. This procedure gives the same answer (|13[) . 
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The matter counterterms do not contribute to the finite part of the action for d > 2. In order 
to evaluate the thermodynamic potentials, the Euclidean "time" direction r is compactified on a 
circle with period At = 4n / (d+2)rQ = 1/((3T). The grand canonical potential f2 = —pV = — TlnZ 
is then given by 

n(T^) = -cVT^ fpj 2 , (14) 

where the (T, /^-independent constant c = (Av )7r d+1 2( 3d ~ 2 ^ 2 (d+2)~( d+2 ^ /Gd+3 determines the 
normalization of all thermodynamic functions, and therefore effectively counts the degrees of free- 
dom in the dual field theory. 4 All other thermodynamic functions can be easily determined as 
n = dp/dfi, s = dp/dT, e = —p + Ts + /in, and one finds 2e = dp, as is expected in a scale- 
invariant theory with z = 2. The dependence on temperature and chemical potential is precisely 
of the expected form (|3|), with the scaling function g(x) = 4p \x\~( d+2 ^ 2 . Note that negative /i is 
required so that the density n = dp/dfi is positive. The total entropy 

S = ~§ = ~n = -(d+2)S E = ^-t^PAvV (15) 
is equal to the Bekenstein-Hawking entropy ([7]), a consistency check for thermodynamics. 



V. VISCOSITY 

In a translation-invariant theory with a hydrodynamic regime at long distances, shear viscosity 
can be evaluated using the Kubo formula 

v = _ i im I Im Gf 2jl2 (w, k=0) , (16) 

where G^ 2 12 i s the retarded function of the component T x \ x i of the energy-momentum tensor. The 
relevant component of the metric perturbation in the bulk tp = h 2 can be Fourier transformed 
along the u and x l directions, ip(u, u,x, r) = fduid d k/(2ir) d+1 (p UJ: ] i (r)e~ , ' UJU+ ' lk ' :x ~ (note that u is 
identified with the time in the dual field theory), and only k=0 is needed to find the viscosity. The 
perturbation is ^-independent, as required by the fact that the energy-momentum tensor has zero 



particle number. The fie 
minimal massless scalar 



d tp decouples from other perturbations, and satisfies the equation of the 



la ] . In the black hole background ([5]) , one finds 



^» + \ { X 2 -r d + 2 ) Vjr) + { r d + 2 -r d / 2 ) 2 ^ ir) = °" ^ 

The equation has a regular singular point at the horizon, with exponents a± = ±4^- Evaluation 
of the retarded function requires that one keeps only the wave going into the black hole fl7| . 
which fixes the solution as <£u{ r ) = (r—ro) a ~ F(r), where F(r) is regular at the horizon r = tq. In 



4 In the analogous gravitational description of relativistic field theories, the normalization coefficient c is uniquely 
related to the central charge of the corresponding CFT [Ti| . 
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the hydro-dynamic limit u — > 0, the equation for F(r) can be solved as a power expansion in uj. To 
linear order, one finds 

F(r) = F - 7 ^ln(^) + O^), (18) 



(d+2) r \r-r 

where Fq is an integration constant which determines the overall normalization of the perturbation. 
It is proportional to the boundary value of the field (po(uj) = ^w(^a)> where r\^oo is the large- 
r cutoff. The two-point function is determined from the regularized on-shell action f)12[) . The 
Gibbons-Hawking term ensures that the action for ip has the form (up to contact boundary terms) 



(W7rG d+3 )S = ~ Jd d+3 x^g-g^d^dv<p 



' ftL: d '' k r d +*f(r)^Ur)^(r)Av\. (19) 



2 J 2vr (2vr) d " v ' r ~ v ' rajv ' ' lr=r A 

Taking the second variation of the action with respect to <po(ui) [ltil . one finds the viscosity 
from the Kubo formula (I16p . using the solution (1181) . 

lb7rG d+3 

Comparing with the expression for the black hole entropy ([7]) , the ratio of shear viscosity to entropy 
density s = S/V is equal to (restoring K) 

VI. DISCUSSION 

Our model ^ can be viewed as an extension of the d = 2 model of Ref. 0] to arbitrary dimension. 
Such a generalization is non-trivial because it involves both guessing the action and finding the 
solution to the equations of motion. We have not attempted to find a string theory embedding of the 
action (jj]) , and we do not know whether black hole solutions ([5]) uplift to ten- or eleven-dimensional 
supergravity for d > 2. Nevertheless, we feel that the holographic model ([!]) is interesting in its 
own right, in particular because it provides a window to a novel class of scale invariant 3 + 1 
dimensional non-relativistic models with strong quantum fluctuations. One way to see that the 
quantum fluctuations in the dual field theory are strong is to note that the viscosity is small 
and saturates the conjectured quantum-mechanical bound of Ref. [3]. This means that the dual 
quantum system can not be described as a gas of weakly interacting quasiparticles because the 
latter would imply a parametrically large viscosity to entropy ratio. Another way to see that 
interactions are important is to note that the thermodynamic scaling function g(x) defined by 
Eq. ([3]) has different form in a non-interacting gas, and in the present holographic model. In a 
non- interacting gas of fermions or bosons in 3 + 1 dimensions, the scaling function is proportional 
to g(x) oc Jq 00 dun 3 / 2 (e u ~ x ± 1) _1 18], while in the holographic model g(x) oc |x| _5//2 . 



The energy per particle E/N = e/n in the holographic model is 

E 3, 
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(in arbitrary dimension, E/N = \ jj\d/ {d+2) ) at any temperature, including T — > 0. Expressed in 
terms of the number density, the energy per particle is proportional to n -2 / 7 ! 110 / 7 (in arbitrary 
dimension, ~ n - 2 /(rf+4)y2(d+2)/(d+4)^ and ag a regult ^yjy _^ o as T ^ at fixed density. For 

unitary fermions at T = 0, it is conventional to parametrize the energy per particle as E/N = ggp £, 
where e F ~ n 2 / 3 is the Fermi energy in a gas of non- interacting fermions. Lattice simulations 19] 
suggest that the thermodynamic scaling function for the black holes discussed in this paper differs 
from that of unitary fermions, in particular giving £ about 0.4 [20]. In the holographic model, 
on the other hand, the temperature dependence of E/N suggests that £ is zero in any dimension. 
However, it is possible that at lower temperatures in our model there exist different backgrounds 
with the same asymptotics and different thermodynamic properties which are closer to those of 
unitary fermions. This would be natural as the background ([5]) is invariant under f-translations. 
Therefore it does not break particle number symmetry and should not correspond to a superfluid 
state in the boundary theory. It would be interesting to find physical systems with the same 
thermodynamic scaling function as in the holographic model. 

In this paper, we only considered a holographic model corresponding to a quantum field theory 
with the dynamical exponent z = 2. It would be interesting to find black holes corresponding to 
non-relativistic systems with z ^ 2 in arbitrary d. Another question to understand is whether 
black holes (0) can be embedded in string theory. Finally, we see that for the non-relativistic field 
theories studied here, the ratio of shear viscosity to entropy density has a universal value of h/4n 
in any dimension d 2, despite the fact that the black hole backgrounds ([5]) do not satisfy the 



assumptions of universality theorems of 15|, [2l|, |22| . This suggests that the universality of shear 



viscosity can be proven for a wider class of black hole spacetimes, and therefore for a wider class 
of field theories. 
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